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SUMMARY 

9 IS*® 

Hemispherical emittance, both total and normal, were calculated from normal 
spectral-emittance data obtained from other sources. The metals evaluated were clean, 
polished tungsten, molybdenum, and tantalum, each of which exhibits spectral emittances 
that vary considerably with temperature and wavelength. 

Net radiant heat flow between two parallel, infinite plates was computed by summing 
the monochromatic energy exchange. The evaluation was made for all nine possible 
combinations obtained by interchanging metals on the two surfaces. The results are 
graphically presented as a function of the temperatures of the two surfaces. Equations 
of the form 



were fitted to each of the nine sets of heat-flux calculations, where q is the heat- 
transfer rate, and T^ and Tg are the temperatures of the hotter and cooler surfaces, 
respectively. Values of the constants, a, b, and c are presented along with contour 
plots showing the temperature regions in which the equations are accurate to 7 and 11 
percent, respectively. These equations are shown to be more simple and considerably 
more accurate than conventional gray- body equations in representing machine- calculated 
nongray-body heat flux. 

The computational procedures for obtaining the heat flux and the constants in the 
arbitrary equations are described. 

INTRODUCTION 




Advances in high-temperature technology during the last decade have resulted in 
many applications requiring refractory metals. These metals are finding uses in their 
pure form, in alloys, and in certain compounds such as carbides that exhibit many metal- 
lic characteristics. The spectral emittance of these materials is highly sensitive to 


wavelength. For a spectrally selective material, the total emittance and total absorp- v 
tance are not equal; hence the conventional gray-body heat-exchange equations based on 
the interchangeability of emittance and absorptance are invalid. The proper procedure 
is to sum the monochromatic energy exchange over all wavelengths of heat radiation. 

This summation method is used herein to evaluate the net heat exchange between two 
parallel, infinite surfaces of clean, polished tungsten, molybdenum, and tantalum. The 
evaluation was made with a digital computer for all nine possible combinations of hot and 
cold surfaces. Each surface was at a uniform temperature that was varied from about 
the melting point to almost 500° K. The results are presented in graphical form and in 
the form of fitted equations. 

Hemispherical spectral and hemispherical total emittances, computed from the nor- 
mal spectral -emittance data of reference 1, are also included. These data, which are 
for very pure, single crystal metals, are compared to results of other investigators. 


ANALYSIS 

The equation and a computational procedure for evaluating the nongray-body heat 
flux between the two parallel, infinite surfaces of opaque material are described in de- 
tail in references 2 and 3. 

The basic method, most frequently called the summation method, incorporates the 
following assumptions or conditions: 

(1) Each surface is at a constant temperature. 

(2) Each surface is a diffuse radiator. 

(3) At a given temperature and wavelength the hemispherical spectral absorptance 
equals the hemispherical spectral emittance. 

(4) A vacuum exists between the two plates. 

The nongray-body net heat flux q n ^ in watts per square centimeter is given by 


q 


ng 




( 1 ) 


(symbols are defined in appendix A. ) In the finite difference formed required by the 
machine calculators, the spectral emittance of wavelengths greater than 20 microns was 
assumed to be constant and the energy at wavelengths less than 0. 2 micron was negli- 
gible. Equation (1) then becomes 
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Ratio of hemispherical to normal emittance 


A=20 



4 

The flux per square centimeter emitted by a surface is ae,pT , where €rp, the 
hemispherical total emittance at a particular temperature, is given by 
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(3) 


In finite difference form, equation (3) becomes 



(4) 


The size of the increment AA was chosen to be constant at 0. 05 micron. 

The hemispherical spectral-emittance data to be used in the foregoing equations were 
computed from experimental measurements of normal spectral emittance (ref. 1). These 

normal emittances are shown in table I (p. 4). 
For tungsten and tantalum, metal tempera- 
tures were 1370°, 1920°, and 2480° K. Mo- 
lybdenum data were obtained only at the two 
lower temperature values. The procedure of 
preparing these data for heat- flux calculations 
is as follows: 

First, the tabulated normal data were con- 
verted to hemispherical spectral emittancps 
by using figure 1, which was obtained by the 
E. Schmidt and Eckert method detailed in ref- 
erence 4. Next, the results were plotted in 

Figure 1. - Ratio of hemispherical to normal emittance of 
metallic surface. 
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TABLE I. - NORMAL SPECTRAL- EMITTANCE DATA OF 


REFERENCE 1 


Wavelength, 



Surface temperature, 

T, °K 



A, 

M 

1370 

1920 

2480 

1370 

1920 

2480 

1370 

1920 


Normal spectral emittance, T n 


Tungsten 

Tantalum 

Molybdenum 

0.42 

0.478 

0.463 

0.456 

0. 520 

0.511 

0. 502 

0.406 

0.388 

. 54 

.462 

.446 

.438 

.494 

.471 

.459 

.392 

.375 

.68 

.445 

.428 

.416 

.426 

.412 

.407 

.374 

.353 

1.00 

.377 

.365 

.356 

.255 

.291 

.312 

.301 

.289 

1.28 

.314 

.309 

.309 

. 196 

.238 

.267 

.233 

.240 

1.60 

.249 

.257 

.266 

. 167 

.205 

.236 

. 179 

. 197 

2. 00 

. 181 

.204 

.224 

. 149 

. 183 

.210 

. 141 

. 167 

2. 50 

. 135 

. 164 

. 191 

. 136 

. 164 

. 189 

. 117 

. 145 

3.00 

. Ill 

. 140 

. 169 

. 126 

. 153 

. 178 

. 103 

. 130 

4.00 

.0879 

. 116 

. 143 

. 114 

. 137 

. 158 

.0895 

. 113 

5. 00 

.0763 

. 102 

. 127 

. 106 

. 126 

. 143 

.0793 

. 102 

6.00 

.0693 

.0945 

. 116 

.098 

. 117 

. 133 

.0716 

.0956 

7.00 

.0642 

.0866 

. 109 

.0919 

. 110 

. 127 

.0676 

.0892 

8.00 

.0606 

.0821 

. 104 

.0884 

. 104 

. 120 

.0647 

.0840 

9.00 

.0568 

.0779 

. 101 

.0825 

.0989 

. 116 

.0589 

.0798 

10. 00 

.0547 

.0753 

.0969 

.0795 

. 0949 

. 114 

.0577 

.0766 

11.00 

.0526 

.0717 

.0945 

.0745 

.0918 

. 109 

.0561 

.0739 

12.00 

.0475 

. 0669 

.0877 

.0736 

.0870 

. 105 

.0526 

. 0696 

13.00 

.0454 

.0677 

.0890 

.0680 

.0876 

. 103 

.0511 

.0687 

14.00 

.0463 

.0658 

.0847 

.0734 

.0839 

. 102 

. 0529 

.0690 

15.00 

.0458 

.0653 

.0847 

.0535 

.0824 

. 105 

.0442 

.0636 


a manner similar to that shown in figure 2(a), and extrapolated to wavelengths of 0. 2 and 
20 microns. Then at each wavelength to be used in the machine input, the data were 
plotted as emittance against temperature. Families of curves faired through the data 
points determined rp at the values of temperature required by the calculation pro- 
cedure. This technique presented a problem. The normal emittance of molybdenum 
had been evaluated at only two temperatures. This situation was handled by assuming 
molybdenum emittance possessed a nonlinear variation with temperature similar to that 
of tungsten and tantalum. As shown in appendix B, the procedure of preparing the 
e. T data provided room for only a small amount of "engineering error”. 

The machine data input were taken at the values of wavelength and temperature 
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Hemispherical spectral emittance, T 



Wavelength, \ p 
(c) Tantalum. 

Figure Z - Spectra l-emitta nee data. 



0 1000 2000 3000 4000 

Surface temperature, T, °K 

Figure! - Hemispherical total-emittancedata. 


listed in table It (p. 6). For any 
given temperature, a linear interpola- 
lation method was used to obtain 
e-y rp at all intermediate values of 
wavelength. Lagrangian interpolation 
or extrapolation was used to obtain 
^ at the desired temperature. 

RESULTS AND DISCUSSION 

The calculated hemispherical 
spectral -emittance data for tungsten, 
molybdenum, and tantalum are pre- 
sented in table EL These values are 
plotted in figure 2 for three of the 
five tabulated temperatures. Hemi- 
spherical total emittance results for 
the three metals are shown in fig- 
ure 3. 

The metal specimens for which 
the original data were obtained 
(ref. 1) were of high-purity mono- 
crystalline structure. The surfaces 
had been polished either chemically or 
electrolytically to remove damaged 
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TABLE n. - HEMISPHERICAL SPECTRAL-EMITTANCE DATA 
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‘material contained in the surface layers. The emittance measurements had been carried 
out in a high-vacuum environment and in accordance with good spectrographic practices. 
Crystal orientation was found to have little influence on emittance results. For each 
metal, the spectral emittance values were smaller than those of unpolished specimens of 
poly crystalline material (ref. 5). This effect was most noticeable at the longer wave- 
lengths of radiation. 

A comparison of the present emittance data with similar data from other sources 
would be of interest. Recent studies (ref. 6) by Parker and Abbott on clean specular 
surfaces of poly crystalline tungsten and tantalum have yielded e,p data at temperatures 
of 1100° to 3000° K. These data can be compared directly to the data shown in figure 3 
(p. 5). At temperatures above 1400° K, the agreement was excellent; however, as the 
temperature was decreased from 1400° to 1100° K, Parker and Abbott reported progres- 
sively smaller emittances than those presented herein. This effect was approximately 
twice as great for tantalum as for tungsten. Maximum observed discrepancies were 
approximately 8 percent. 

The nongray-body heat-transfer data for the three metals are presented in figure 4. 
Net heat transfer is plotted as a function of T^ - T 2 with T^ as the parameter. The 
heavy solid lines show the region encompassed by the smoothed data from reference 1. 
The lighter solid lines show the calculations based on the temperature extrapolations 
previously described. The dashed lines are the values beyond the temperatures of the 
machine calculations and are simply faired-in curves of heat transfer as a function of 
Tj - Tg. It is well to note that rather large heat-flux errors could accrue from ex- 
trapolation of the e. rn data to the lower temperatures shown in figure 4. 

A , 1 

All data curves of figure 4 exhibit a maximum in the value of q n g. This maximum 
is most noticeable under the combined effect of large T^ and a tantalum cold surface 
(figs 4(c), (f), and (i)). One contributing factor is the extent the absorptance of the 
colder surface decreases with decreasing temperature. At wavelengths near 1 micron 
(where an appreciable amount of blackbody radiation occurs at high temperature levels), 
only tantalum exhibits a decrease in spectral absorptance as the temperature is lowered 
(see fig. 2, p. 5). 

The assumption that the surfaces radiate diffusely can lead to error in the values of 
q n g. An alternate method (ref. 7) of calculating q n g assumes the energy exchange is 
entirely specular and that the angular distribution of energy is defined by the Fresnel 
equations. Exploratory calculations indicate that this alternate method would yield q 
values 3 to 5 percent greater than those presented herein. Since the radiant energy ex- 
change between real surfaces is neither entirely specular nor entirely diffuse, it is 
fortunate that the two limiting situations yield such good agreement in the value of net 
heat flux. 
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(a) Hotter surface, tungsten; colder surface, tungsten, tb) Hotter surface, tungsten; colder surface, molybdenum. (c) Hotter surface, tungsten; colder surface, tantalum. 

Figure 4. - Net radiant heat flux between two infinite, parallel plates. 






Temperature difference, Tj_ - °K 

(b) Hotter surface, tungsten; colder 
surface, tantalum. 

Figure 5. - Comparison of gray- and 
nongray-body methods for computing 
net heat flux between two infinite, 
parallel plates. 


TABLE III. - CONSTANTS OF FORCE-FIT EQUATION (EQ. (6)) 


% - 



Hotter 

surface, 

T i 

Cooler 

surface, 

T 2 

Constant 

Temperature 
of hotter 
surface, 

T i- 

°K 

Equation 
(fig- 6) 

a 

b 

c 

Tungsten 

Tungsten 

0.353 X 10' 15 

5.0 

0 

<3000 

a 


Tungsten 

. 166 X 10" 14 

4.8 

0 

>2000 

b 


Molybdenum 

.314 X 10" 15 

5.0 

0 

<3000 

c 


Molybdenum 

. 149 x 10" 14 

4.8 

0 

>2000 

d 


Tantalum 

.335 X 10" 15 

5.0 

. 17 


e 

Molybdenum 

Tungsten 

0.316 X 10' 15 

5.0 

0 


f 


Molybdenum 

.286 

5.0 

0 


g 


Tantalum 

.308 

5.0 

. 17 


h 

Tantalum 

Tungsten 

0. 328 xl 0~ 15 

5.0 

0 


i 


Molybdenum 

.294 

5.0 

0 


j 


Tantalum 

.326 

5.0 

. 17 


k 
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SIMPLE EQUATIONS FOR NET HEAT EXCHANGE 


Simple expressions for net heat flux ease the mathematical complexity of analytical 
investigations. Two general types of equations were examined; namely, those based on 
conventional gray-body heat flow and those that express the net heat flux as an explicit 
function of and T 2 . 

The conventional gray-body net heat-transfer equation for closely spaced surfaces 
is 



For this equation to be correct, the total emittance of a surface must equal its total ab- 

sorptance. Since the refractory metals do not adhere to this restriction, discrepancies 

between the results of this equation and equation (1) are very large. The correlation can 

be improved by an empirical method (ref. 8, p. 63), which replaces e T in equation (5) 

l 2 

by an effective emittance of the colder surface erp* evaluated at the geometric mean 
temperature 

T* = yTj_ X t 2 

The inaccuracy of this empirical equation is illustrated in figure 5(a), which is typical of 
those cases wherein tungsten or molybdenum was assigned to the colder surface. The 
nongray-body flux is larger than the gray-body flux and, in general, the error increases 
as Tj is decreased and/or as AT is increased. Where tantalum was the colder sur- 
face, typical results are shown in figure 5(b). 

Expressions that would possibly provide better correlations than those demonstrated 
in figure 5 were examined. Of the several expressions investigated, the "force-fit" 
equation 

%-( T J-T^) (6) 

yielded flux values most closely matching the values of equation (1). The constants a, 
b, and c for each set of metal surfaces were arrived at in the following manner: Af- 
ter assigning c a value of zero, a best value for b was determined. This was 


i 


12 


-accomplished by fixing b at each of a number of set values between 4 and 6, and then 
solving for a while ranging the variables and Tg. The best value of b was the one 
that produced a minimum deviation in the value of a. With b known, a best value for 
the constant c was determined in a manner similar to the determination of b. Values 
of c smaller than 0. 1 had little influence and therefore were assigned a value of zero. 
Finally, with the constants b and c established, the constant a was picked such that 
a close fit would be obtained over the widest useable range of temperatures. 

There exists an unlimited number of procedures for selecting the constants. The 
procedure just described is better than some and is easily prepared for use in machine 
calculators. The results are presented in table III (p. 11), and contour plots illustrating 
equation accuracy are shown in figure 6. When constants were being fitted to equa- 
tion (6), difficulties occurred in obtaining a reasonably satisfactory fit for several cases 
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(a) Hotter surface, tungsten. 



Error, 

percent 



] 7 < Error <11 




(b) Hotter surface, molybdenum. 



Temperature of hotter surface, Tj, °K 

Cooler surface: Tungsten Molybdenum Tantalum 

(c) Hotter surface, tantalum. 

Figure 6. - Error incurred with use of force-fit equations from table III. (Error 
is measured in percent deviation of net heat flux obtained by force fit from net 
heat flux obtained by summation. ) 
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in which tungsten was the hotter surface. In each of these situations, two sets of con- 
stants are presented. 

The values of flux given by the force-fit equation were larger than the nongray-body 
values in the low temperature region. No attempt was made to improve the fit because of 
the lack of low-temperature emittance data and hence the possible inaccuracy of the 
nongray-body heat-flux data. 

Comparison of figures 5 and 6 shows that the force-fit equations more accurately 
represent the results of nongray-body heat-transfer calculations than the empirical gray- 
body equations. Also, within the larger shaded regions of figure 6 (where the error is 
less than 7 percent), it is very likely that the q^ calculations are as accurate a mea- 
sure of net heat flux as can be achieved in practice for real surfaces in a parallel-plate 
geometry. 

The results of table III taken as a whole, show that the net heat exchange between 
two parallel, infinite surfaces of typical refractory metals is closely proportional to the 
fifth power of temperature throughout the temperature range investigated. The 
temperature- ratio factor is significant only for those cases in which tantalum is the 
cooler surface. 


RESUME 

Hemispherical emittance, both total and normal, were calculated from normal 
spectral-emittance data obtained from other sources. The metals evaluated were clean, 
polished tungsten, molybdenum, and tantalum. Net radiant heat flow between two closely 
spaced plates was computed by summing the monochromatic energy exchange. This eval- 
uation was made for all nine possible combinations obtained by interchanging metals on 
the two surfaces. 

Simple equations were M force-fit M to the results of the heat-flux calculations. 
Reasonably good fits were established for a wide range of temperatures. The force-fit 
equations show much better agreement with machine -calculated heat fluxes than conven- 
tional gray-body equations, and are also easier to use. 

Lewis Research Center, 

National Aeronautics and Space Administration 
Cleveland, Ohio, March 17, 1965. 
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APPENDIX A 


SYMBOLS 


a, b, c constants 

T Planck's distribution function, 

37 404A“ 5 (e 14 387 / AT . i)" 1 , 
W/(M)(sq cm) 

q heat-transfer rate 

q ff net heat-transfer rate for 

"force-fit" equation, 

W/sq cm 

q^ net heat-transfer rate for gray- 

body conditions, W/sq cm 

q ng net heat-transfer rate for 

nongray-body conditions, 
W/sq cm 


T surface temperature, °K 

T* yTj X t 2 

e T hemispherical total emittance 

rp hemispherical spectral 

emittance 

e > t n normal spectral emittance 

Aj J- y n 

A wavelength, ju 

a Stefan- Boltzmann constant, 

5. 6699x10“ 12 W/(sq cm)(°K 4 ) 

Subscripts: 

1 hotter surface 

2 cooler surface 


15 



APPENDIX B 


EXTRAPOLATION OF EMITTANCE DATA 

Equations (1) and (2) require the emittance be independent of emission angle, that is, 
each surface emits according to the cosine law of Lambert. Therefore the normal emit- 
tance data of table I had to be converted to hemispherical values, ,p. For diffuse 
radiators, a one to one relation exists. For polished metal surfaces, the hemispherical 
emittance can be obtained by means of figure 1. This conversion curve is based on equa- 
tions (4-79), (4-84), and (4-85) of reference 4 relating electrical resistivity to emit- 
tance. The equations of little utility at the shorter wavelengths, appear satisfactory at 
wavelengths greater than 4 microns for smooth tungsten surfaces (ref. 2 and fig. 22, 
ch. 13 of ref. 9). A comparison of figures 1 and 2 shows the ratio q/ e A T n to 
largest in the region of greatest reliability of the equation. Hence, for polished sur- 
faces, the conversion curve provides a fairly reliable means of obtaining heat-transfer 
results based on normal spectral-emittance data. 

The aforementioned equations of reference 4 were used as guidelines to extrapolate 
the hemispherical-emittance data to wavelengths beyond 15 microns and to conduct tem- 
perature extrapolation of the emittance data at wavelengths greater than 4 microns. 

Since very little energy is transferred at wavelengths less than 0. 42 micron, a sim- 
ple extrapolation, based on data of Larrabee (ref. 10) was used to extend the hemi- 
spherical emittances to a wavelength of 0. 2 micron. Temperature extrapolation at short 
wavelengths was based on observations (ref. 9 and elsewhere) that for temperatures at 
least as small as 300° K, the emittance data for refractory metals pass through an in- 
flection point (fig. 2). 
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